Detailed Atmosphere Model Fits to Disk-Dominated ULX Spectra 
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ABSTRACT 



We have chosen 6 Ultra-Luminous X-ray sources from the XMM-Newton 
'. archive whose spectra have high signal-to-noise and can be fitted solely with 

a disk model without requiring any power-law component. To estimate system- 
atic errors in the inferred parameters, we fit every spectrum to two different disk 
models, one based on local blackbody emission (KERRBB) and one based on 
q ! detailed atmosphere modelling (BHSPEC). Both incorporate full general rela- 

£j ! tivistic treatment of the disk surface brightness profile, photon Doppler shifts, 

and photon trajectories. We found in every case that they give almost identical 
fits and similar acceptable parameters. The best-fit value of the most interesting 
parameter, the mass of the central object, is between 23 and 73 M in 5 of the 
6 examples. In every case, the best-fit inclination angle and mass are correlated, 
in the sense that large mass corresponds to high inclination. Even after allowing 
for this degeneracy, we find that, with > 99.9% formal statistical confidence, 3 of 
the 6 objects have mass > 25 M ; for the other 3, these data are consistent with 
a wide range of masses. A mass greater than several hundred M is unlikely for 
the 3 best-constrained objects. These fits also suggest comparatively rapid black 
,£h ! hole spin in the 3 objects whose masses are relatively well-determined, but our 

estimate of the spin is subject to significant systematic error having to do with 
uncertainty in the underlying surface brightness profile. 

Subject headings: black hole physics - accretion disk - IMBH - X-ray 



Introduction 



Nearly two decades since their first unveiling in the X-ray sky (IFabbiano &: Trinchieri 



19871 ). Ultra-Luminous X-ray sources (ULXs) still remain puzzles. The X-ray luminosities 
of these off-galaxy-center point sources fall in the range 10 39 - 10 41 erg s _1 in the 2-10 keV 
band, at least an order of magnitude larger than those of Galactic X-ray Binaries (GXB) and 
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several orders of magnitude smaller than those of Active Galactic Nuclei (AGN). If ULXs 
are accreting black hole systems that observe the Eddington limit (L# = 1.3 x 10 38 erg s _1 ), 
the central black holes must have masses from a few tens to a few thousands of solar masses, 



Colbert & Mushotzkv 


1999; 


Miller & Hamilton 


2002; 


Miller et al. 



20031 ). But it is also possible that ULXs' masses are smaller than the Eddington limit 
suggests. They might be st ellar mass black holes whose radiation is beamed in the observed 



direc t ion (IKing et al. 



2001) or which are accreting at super-Eddington rates (jBegelmanll2002 



20061 ; IStern &; Poutanenll2006l ). The only way to distinguish between these scenarios is to 
measure the central masses of ULXs as accurately as possible. 

The most reliable way to measure the masses of the stars in a binary is to observe 
their orbital motions. Unfortunately, this method does not work in the context of ULXs. 
Because there are so far no clearly identified spectral lines in ULXs' spectra, the com- 
mon technique used in spectroscopic binaries cannot be applied. Temporal variability is 
another approac h. After compa ri ng th e Fourier power spectra of GXBs and AGNs, some 



researchers (e.g. iMcHardy et al.l (120061 )) have concluded that the characteristic timescales 
of accreting black holes are proportional to their central objects' masses. If ULXs are ex- 
tragalactic analogs to the GXB systems in our galaxy, similar correlations between quasi- 
periodic oscillation (QPO) frequencies and mass may be used to estimate ULX masses. 
These correlations have been est ablished successfully for several ULXs, for ex ample, M82 XI 
dstrohmaver k Mushotzkvll2003h and NGC 5408 XI Jstrohmaver et alil2007h . The discovery 
of the QPO in NGC 5408 XI led to an estimated ULX mass of 1500 - 3500 M . However, 
these works treated a limited number of ULXs and it is not guaranteed that QPOs can be 
found in a large sample of ULXs. Even if there is a QPO identified in the power spectrum 
of a ULX, the mass inferred through this method is hi ghly dependent on the "calibration 
standard" and may be quite uncertain (see, for example IStrohmayer et al.l 120071 ) . 



Another method for inferring the central mass is spectral fitting (jVierdayanti et al.ll2006 



Winter et all2006l . 120071 ). To fit a ULX spectrum, various models and/or model combinations 
can be exploited. In practice, ULXs' spe ctra are usually modelled by a combination of a 
thermal disk component and a power-law (IMiller et al.ll2003l . 120041 ; iRoberts et al.l 120041 )). In 
this combination, the thermal disk component covers the soft X-ray (0.3 - 1 or 2 keV) range 
of the spectrum while the power-law takes care of the hard tail (1 - 10 keV and above). 
Because the spectrum of soft X-ray radiation from the accretion disk is dependent on the 
central mass, but there is no clear relation between the hard tail of the spectrum and the 
central mass, we define our sample to include exclusively those ULXs with only thermal 
disk components in their spectra. In this way, we can focus purely on the component of the 
system — the thermal accretion disk — with the tightest connection to the central object's 
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mass, as well as the other key accretion parameters: spin, accretion rate and inclination 
angle. 

There are many thermal disk models that have been used in this way In Xspecll, 
for instance, there are DISKBB, DISKPN, GRAD, KERRBB. Thes e models share more 
or less the same core — the multi-color disk (M CD) approximation (IMitsuda et al.l Il984j ; 



Makishima et al.lll986l : IShimura fc Takaharalll995l ). In this approximation, one assumes that 
the disk emission is a sum of a series of local blackbody spectra with effective temperatures 
defined by the standard disk model but modified by a hardening factor. The version we 
employ here is KERRBB, which also incorporates genera l relativistic effect s. Even though 
the MCD approximation works well in the case of GXB (IDavis et al.ll2005l ). in this paper, 
we also fit the data to a model in which a detailed stel lar atmosphere is computed for each 
disk annulus. Our reason to do so is that earlier work (IHui et al.ll2005l ) has shown that the 
spectral shape, especially around the peak energy, can be altered by atomic features when 
the central mass is larger than typical stellar masses (> 100 M ). This model (which also 
fully includes general relativistic effects ) has b een compiled to an easy-to-use table model 
in Xspecll (BHSPEC: IDavis Hubenyl (120061 )): we modified it and extended its parameter 
space to cover a wider range of central masses (see § 4 for details). 

In § 2 we describe the sample selection criteria and in § 3 the data reduction procedure. 
The disk model BHSPEC is briefly discussed, and specific fitting results for each object are 
presented in § 4. For comparison, previous work on the objects in our sample is briefly 
described in § 5. We discuss our results in § 6 and draw conclusions in 5 7. 



2. Sample selection — disk dominated ULXs 



Our sample was drawn from the XMM-Newton public data archive . We began with 
the object lists compiled by Winter et. al. (see Table 6 in IWinter et al.l (120061 ) and Table 



1 in 



Winter et all (120071 1) 



For their sample, Winter et al. chose objects observed for at 
least 10 ks and with distance less than 8 Mpc. These criteria guaranteed a minimum of 
400 counts for objects with Lx > 2 x 10 38 erg s _1 . To focus on ULX candidates, they also 
removed all objects with spectra distinctly different (e.g., resembling supernova remnants) 
and all objects located at the centers of their host galaxies (to eliminate AGN). We chose 
those ULX candidates with > 1000 counts in order to guarantee good signal-to-noise in the 
reduced spectra. Then, in order to retain as much as possible of the soft X-ray emission 
from the accretion disk, we excluded those objects with neutral hydrogen column density 
in the line of sight larger than 10 22 cm" 2 as dete r mined by Winter et al. (see Table 4 in 



Winter et all (120061 ) and Table 2 in IWinter et al.l (l2007h ). except for Circinus XI and X2 
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which have slightly higher values. There are 23 objects meeting these criteria. 

We then used the latest version (7.0) of the XMM-Newton Science Analysis System 
(SAS) and calibration (up to date as of Jan. 2007) to reprocess the data and obtain the 
event files. After applying proper filters (see § 3 for details), source and background spectra 
for each of these objects were created and fitted by the "disk+powlaw" combination in 
Xspecll. Specifically, we used a photoelectric absorption model (phabs) to account for the 
opacity in the interstellar medium and the thermal disk model BHSPEC (a detailed disk 
model discussed later in this paper) plus a power-law component. To decide which model 
component (disk or powlaw) dominates the spectral output, we computed the integrated 
counts from both components and defined the spectrum to be dominated by one or the 
other when one component contributes more than 70% of the total counts. There are 11 
powlaw-dominated spectra, 6 disk+powlaw spectra, and only 6 objects (see Tabled]) having 
disk-dominated spectra. In other words, we can use the disk model alone to fit those 6 disk- 
dominated spectra. In these objects, a non-zero power-law contribution at best marginally 
improves the quality of the fit and sometimes even harms it. 



3. Data reduction 

We used the latest version (7.0) of SAS and the up-to-date calibration (Jan. 2007) to 
reprocess the Observation Data Files (ODF). Commands "emchain" (for EPIC-MOS) and 
"epchain" (for EPIC-PN) in SAS 7.0 were used to get the new processed pipeline products 
(PPS) which include the photon event files. Standard data filter procedures were followed 
according to the instructions in the XMM-Newton ABC guide. Good events were required 
to satisfy the following conditions: 1) for the MOS (both MOS1 and MOS2) detectors, the 
event pattern is in the to 12 range (single, double, triple, and quadruple pixel events) and 
the pulse height is in the range of 0.2 - 12 keV; 2) for the PN detector, the event pattern 
is in the to 4 range (single and double pixel events) and the pulse height is in the range 
of 0.2 - 15 keV. The #XMMEA_EM (for the MOS) filter and #XMMEA_EP (for the PN) 
filter were also applied together with a "FLAG==0" filter to kick out any bad or close-to- 
edge events captured by the CCDs. The light curve for each observation was produced with 
the SAS command 'evselect" to decide a proper time filtering threshold in order to reduce 
the influence from high rate flaring. For the MOS detectors, the good time intervals were 
selected by setting the rate to less than 5 cts s _1 (sometimes a little bit higher); for the PN 
detector, the threshold varies from 25 to 60 cts s _1 . 

After obtaining the filtered event files, source and background spectra were extracted 
by using the SAS procedure "especget". For the sources, circular regions with radii of 20 
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arcseconds (±10 according to the size of the sources and closeness to other sources or detector 
edges) were applied in the spectra extraction. For the background, a nearby circular region 
close to but without overlap on the source's region was used, with radius double the size 
of its corresponding source. SAS procedures "rmfgen" and "arfgen" were used to generate 
the response matrix files (RMF) and ancillary response files (ARF). Finally, we used the 
command in HEASOFT 6.1 — "grppha" — to regroup the spectra so that there are at least 
20 photon counts in every spectral bin. 



4. Model description and fitting results 
4.1. BHSPEC vs KERRBB 



We fitted the selected spectra with both BHSPEC (a detailed atmos 



phere calculation) 



and KERRBB (a MCD model). For details of KERRBB, see lLi et all (120051 ) and the Xspecll 



manual; the only significant choice we made was to set th e dilution factor to 1.7. B HSPEC 
is based on an atmosphere model of a standard a— disk (IShakura fe Sunyaevl Il973l ) where 
each ring is treated in hydrostatic equilibrium and energy balance. The stress parameter a 
is chosen to be 0.01 so that all disk annuli are optically thick. The radiation intensity is 
computed as a function of both frequency and position by means of a full radiation transfer 
solution. Continuum opacities due to free elect rons, H, He, C, N, O, Ne, Mg, Si, S, Ar, 
Ca, Fe and Ni, with the abundances given by lAnders fe Grevessd (119891 ) are calculated 
on the basis of explicit ionization balance calculations and statistical equilibrium of the 
most populated internal states. Photon-electron energy exchange by Comptonization is also 
included in the transfer solution. General relativity effects in photon propagation are taken 
into account: radiation is boosted and beam ed, and its trajectories "bent" in t he general 



relatiy istic potential before it reaches infinity ( lAgollll997l ). For further details, see lHui et al. 

(boosh . 



In the work of iDavis fe Hubenyl (120061 ). disks with many black hole masses, spins, ac- 
cretion rates, and inclination angles were modelled, and the spectra obtained were compiled 
into the Xspecll table model BHSPEC. For this project, we extended the parameter space 
considered by Davis and Hubeny to include disks around black holes with masses in the range 
of 10 - 10 4 M . Because the internal consistency of disk models becomes dubious when the 
luminosity approaches or exceeds the Eddington luminosity, we included no models with 
L/Le > 1. When we speak of lower bounds on the mass based on BHSPEC-fitting, we 
therefore always mean that this bound is subject to the constraint of being consistent with a 
disk model, i.e., keeping the luminosity sub-Eddington. It should be borne in mind, however, 
that some significant flexibility nonetheless remains at the low-mass end of parameter space 
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because altering the inclination angle or the black hole rotation rate could keep the luminos- 
ity in our direction at the level required by the measurements even while the luminosity in 
other directions is diminished in order to satisfy the sub-Eddington constraint. We did not 
so restrict KERRBB, so any super-Eddington fits it produces must be regarded as extremely 
suspect because the underlying model would then be inappropriate. For technical reasons, 
we limited the range of spins considered by BHSPEC to < a/M < 0.997 and by KERRBB 
to -0.999 < a/M < 0.999. 



4.2. Fitting results 

The fitting results for the objects in our sample are shown in Table [2j Each object is 
given two rows. The first row, with the object's name, describes the fit with BHSPEC; the 
second row shows the result from KERRBB. In BHSPEC, the normalization is fixed at a 
value corresponding to the best estimate for that object's distance (tabulated in Tabled]). All 
other parameters, i.e., the equivalent hydrogen column density of photoelectric absorption, 
the black hole mass and spin, the accretion rate, and the disk viewing angle are set free 
to vary. In KERRBB, the parameter "distance" is fixed to the same value as in BHSPEC 
for each object and the normalization is set to 1.0. KERRBB also offers several other 
parameters. We freeze "eta" to zero (the zero torque inner boundary condition), "hd" to 1.7 
(common diluted black body spectrum), "rflag" to —1.0 (self- irradiation is off), and "lflag" 
to 1.0 (limb-darkening is on). 

The best-fit values and their uncertainties are shown in Table El All are statistically 
acceptable, and in each case the BHSPEC least x 2 is slightly smaller than the KERRBB 
value. All one-parameter uncertainties were computed with Ax 2 = 2.706, equivalent to 90% 
confidence for a single parameter. For the BHSPEC fits, besides the best values, we also 
calculated the ratio of the observed luminosity (E = 0.3 - 10 keV, except for Circinus X2 
MOS1, in which case the integration is limited in E = 0.55 - 10 keV) to the computed total 
intrinsic luminosity (E — - oo keV) emitted in our direction; these are shown in the last 
column of Table [2j The fact that all these ratios are close to unity demonstrates that the 
XMM-Newton energy band contains nearly all the light our models suggest is being radiated. 

In the case of KERRBB, the spectral fitting gives the accretion rate in the form of M 
instead of the dimensionless L/Lg. We used the following formula to translate between those 
two definitions, assuming that the radiative efficiency, rj, is the standard function of a/M 



f lNovikov fc Thorndll973f ). 

L rjMc 2 



L E 1.3xl0 38 M' 
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Here c is the speed of light, M is the accretion rate in units of gm s _1 and M is the accretor's 
mass in units of M . Because the total luminosity of an object is fixed when the distance 
is known, the ratio L/Le is actually inversely proportional to the object's mass. For this 
reason, only the uncertainty in the estimated mass is taken into account in the quoted 
uncertainty for L/Le- 

In the following subsections, we discuss each object individually. The spectra and best- 
fit models (BHSPEC only) are shown in the multipanel Figure [TJ The KERRBB fits are so 
nearly identical to the BHSPEC fits that figures displaying them would be visually indistin- 
guishable from each other. 

4.2.1. M81 XI 

We used a 120 ks XMM-Newton observation (0200980101) of Holmberg IX (made in 
September 2004) for M81 XI, which appears in the same field. We extracted spectra from 
the event files generated for the MOS1 and MOS2 cameras and fitted them with BHSPEC 
and KERRBB. The mass of the black hole is inferred to be 67 - 85 M Q from BHSPEC, while 
KERRBB gives a partially overlapping, but rather wider possible range, 33 - 74 M Q . For the 
accretion rate, both BHSPEC and KERRBB indicate values close to the Eddington limit, 
which makes the applicability of both models suspect because both assume thin disks. Both 
models also give similar values for the hydrogen column density (1.5 - 2.0 xlO 21 cm -2 ), the 
black hole spin (close to the maximum Kerr value) and inclination angle (50° - 70°). 

4.2.2. M101 X2 

M101 was observed (0104260101) in June 2002 for a duration of 43 ks. We extracted 
spectra for the MOS1, MOS2 and PN cameras. The BHSPEC fit suggests that the black 
hole in M101 X2 has a mass 30 - 178 M and the spin of the black hole is close to the 
maximum Kerr value. The accretion rate may be anywhere from the Eddington limit down 
to a moderate level {L/Le — 0.3). The viewing angle is inferred to be 42° - 78°. Using 
KERRBB indicates a lower hydrogen column density and larger uncertainty ranges for the 
mass and accretion rate even though the best-fit values are not far from those found using 
BHSPEC. The spin and the viewing angle in the KERRBB fit are totally unbounded. 
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4-2.3. NGC 253 XI, X3 and X4 

The observation (0152020101) of NGC 253 was made in June 2003 with a duration of 
140 ks. The spectra were drawn from the event files of the MOS1, MOS2 and PN cameras. 

For NGC 253 XI, both BHSPEC and KERRBB indicate similar hydrogen column den- 
sities. For the black hole mass, BHSPEC infers a 90% confidence range of 50 - 81 M , while 
KERRBB gives a very similar, but slightly wider range, 54 - 95 M Q . In both models, the 
spin is quite large and very close to the maximum Kerr value. The Eddington-normalized 
luminosity inferred by BHSPEC is 0.35 < L/L E < 0.56, while KERRBB suggests a lower 
range, 0.18 < L/L E < 0.31. The viewing angle is well constrained by both models in the 
range 64° - 83°. 

In the case of NGC 253 X3, the uncertainty in the black hole mass derived from both 
BHSPEC and KERRBB is somewhat greater, from 25 to 80 M for BHSPEC, from 32 to 
176 M Q for KERRBB. In terms of L/L E , the range preferred by BHSPEC (L/L E from 
0.13 to 0.39) overlaps that of KERRBB (L/L E from 0.03 to 0.18), but for the most part 
suggests higher values. Both models give their best fits when the black hole spins rapidly. 
The viewing angle is large, estimated to be between 44° and 85°. 

The inferred black hole mass from BHSPEC for NGC 253 X4 lies in the range 10 - 
91 M Q ; KERRBB gives a similarly large uncertainty, from 5-77 M . For the normalized 
luminosity, the two models have different best-fit values (L/L E = 0.14 for BHSPEC and 0.05 
for KERRBB), but the ranges permitted within the errors overlap to a considerable extent: 
0.06 - 0.24 (BHSPEC) and 0.02 - 0.35 (KERRBB). BHSPEC suggests rapid black hole spin 
with a large viewing angle (43° - 88°), while KERRBB puts no constraint on either of these 
two parameters. 



4.2.4. Circmus X2 

Circinus X2 was observed (0111240101) in August 2001 for a duration of 110 ks. We 
extracted its spectra for both the MOS1 and PN cameras. Both BHSPEC and KERRBB 
failed to constrain almost all parameters except the hydrogen column density (5.0 - 5.9 xlO 21 
cm -2 ). Unlike the other objects in our sample, the best-fit values of the mass inferred by 
BHSPEC and KERRBB (i.e., 340 M vs. 13 M ) are quite different, while the best-fit values 
of cos 2 are also very different (0.07 vs 1.00). This apparent conflict will be addressed later 
in the discussion when we investigate the correlation between inferred mass and inclination 
angle. 
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5. Comparison to previous work 



In a tho rough study o f M81 XI (there named M81 X6) based on a May 2000 observation 
by Chandra, ISwartz et al.l (120031 ) (hereafter, S2003) showed that the best-fit model for the 
reduced spectrum is an absorbed disk blackbody. This disk domination in the spectrum is 
confirmed in our study. Swartz et. al. found the hydrogen column density (in units of 10 21 
cm" 1 ) to be 2.17 ± 0.10, which is close to our fitting results, Oll^ 2 for BHSPEC and 
1.611qii f° r KERRBB. The black hole mass they estimated was 18 M , but they did so 



assuming the black hole is non-rotating. According to our fitting results (for both BHSPEC 
and KERRBB), the spin of the black hole in M81 XI has close to the maximum Kerr 
value. These results are compatible with our estimate that the mass is 33 - 85 M because 
the (Boyer-Lindquist) radial coordinate of the innermost stable circular orbit (ISCO) for a 
maximal Kerr BH is ~ 4 that of a Schwarzschild BH of the same mass. 



In S2003 and a previous ASCA study of M81 XI (IMakishima et al.ll2000l ). the total 
luminosities were reported to be equal to or even exceeding the Eddington limits of the 
inferred black hole masses. To avoid this violation of conventional expectations, a bigger 
black hole mass (to give a higher Eddington limit) with a larger spin (to make a smaller 
ISCO) was suggested by some authors. As shown by our results, however, if one is constrained 
by fitting the spectrum, there is a limit to how much the mass and spin can be increased, 
and neither our atmosphere-based spectral model nor the multi-color disk model is consistent 
with substantially sub-Eddington behavior. 



Jenkins et al.l (120041 ) studied the same XMM-Newton observation of the galaxy M101 
as we did, but named M101 X2 M101 XMM-1. They fitted the spectrum with a single 
disk component (DISKBB in Xspec) with satisfactory statistics and obtained an inner disk 
temperature (Tj n = 1.33 keV). Due to the limited predictive power of DISKBB, no further 
conclusions were reached in their work. 



6. Discussion 

Before dealing with our principal concern, the inferred masses of these objects, we 
first discuss a more technical point: the surprising lack of difference between BHSPEC 
and KERRBB when fitting the spectra. As described in § 14.11 BHSPEC is a much more 
sophisticated model than KERRBB and one might expect that the atomic features it can 
predict might lead to interestingly different fitting results. The fundamental reason why 
this does not happen is that, at the color temperatures of the objects in our sample, even 
Fe is fully-stripped. It is worth elaborating briefly on this point because this constraint 
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enters in a somewhat indirect fashion. To zeroth order, the temperature of the disk scales as 
L l l A M~ l l 2 . If there were no relativistic Doppler shifts, the color temperature and luminosity 
would then suffice to determine the mass. However, disks around black holes do offer large 
Doppler factors, opening a wider parameter space in which to search for acceptable models. 
For these objects, we find that in part of this parameter space, the mass is high enough 
(> 100 Mq) that the temperature in the disk drops to the point where the unstripped 
fraction of the h eaviest abundant elem e nts can produc e an interesting level of opacity in 



atomic features (IDavis fc Hubenyl 120061 ; iHui et al.l 120051 ). However, the very fact that the 



fluid-frame temperature is this low means that, in order for the spectrum generated to fit 
the data, it must be Doppler boosted, which always implies large inclination angle, and can 
be enhanced by rapid black hole spin. However, the non-uniformity of the boost around the 
disk surface also entails strong smearing of sharp features. The result is that atomic features 
are never apparent in any of our fits, even in the extreme case of Circinus X2, in which 
the best-fit mass is 340 M . Moreover, because the unstripped ion fractions are so small 
at thermodynamic temperatures close to the observed color temperatures, this conclusion 
is only very weakly dependent upon elemental abundances. There is however, one possible 
exception to these arguments: in the presence of some atomic opacity, the error entailed in 
KERRBB by assuming a Planckian output spectrum with a fixed dilution factor may be 
larger than in a case in which there is truly zero atomic opacity. 

We now turn to the main goal of this project: our attempt to use spectral fitting of 
thermal disk spectra to constrain the masses of ULXs. In Figure [H we plot the best-fit values 
and uncertainties of the black hole masses for all our objects except Circinus X2. As that 
figure shows, the best-fit values of the black hole masses span a surprisingly small range, 
between 20 and 80 M Q . In fact, the error bars are large enough that the fits are consistent 
with all five having the same mass, ~ 75 ± 10 M Q . 

Moreover, this range appears to be quite differe n t from that seen in Galactic black hole 



binaries. As reviewed by iMcClintock fc Remillard I (120041 ). of all 17 GXBs with measured 



black hole masses, there are none in which the mass is likely to be smaller than 3 M or 
larger than 18 M . 

However, before we can conclude that these objects have a distinctly different mass 
distribution than the Galactic black holes, it is important to look more close ly at the domi 



nant s ource of uncertainty in the mass: the degeneracy between M and cos i. ISun fc Malkan 



(119891 ) pointed out that these two quantities are related almost inversely. In their disk fitting 



analysis of quasars and Seyfert galaxies, they found 

log(M) oc — bcosi, 

where b is a parameter that varies between 0.6 (for a/M = 0) and ~ 1.2 (for a/M = 0.998). 
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The origin of this relationship lies in the fact that higher inclination produces greater blue- 
shifting of light from matter orbiting toward us, thus compensating for the diminution of 
temperature that generally occurs with larger central mass (and therefore radiating area). 
The coefficient is a function of spin because the orbital speed near and outside the ISCO 
increases with more rapid rotation. 

We find a very similar degeneracy in our fits. It is illustrated in the multipanel Figure [3j 
These figures plot the confidence level contours obtained from fitting to the BHSPEC model 
as seen in the M - cos 2 plane embedded in the higher- dimensional \ 2 parameter space. 
The value of \ 2 as a function of M and cos i is found by minimizing y 2 over all the other 
parameters for those two values of mass and inclination. There is indeed a correlation 
between the mass and the inclin ation angle for e very o bject. In addition, two dashed lines in 



the figures show the slope of the lSun &: Malkanl (119891 ) relations, the steeper one for maximal 



spin, the shallower one for no black hole spin. In most cases, as one might expect, their slopes 
bracket the slope of the correlation. 

This degeneracy explains the large contrast between the BHSPEC and KERRBB best- 
fit masses for Circinus X2. In the former case, the best fit has cosi ~ and M ~ 340 M Q , 
while in the latter, cosi ~ 1 and M ~ 13 M . Nonetheless, in terms of the BHSPEC fitting, 
these two very different answers differ statistically by only 1 - 2er. 

With this degeneracy in mind, we now return to evaluating the fitting results for the 
other five objects in our sample. If one assumes a disk model, three of these — M81 XI, 
M101 X2, and NGC 253 XI — require relatively large masses, in all three cases > 25 M Q . 
Smaller masses would lead to very large x 2 , no matter what other parameters (spin, inclina- 
tion angle, intervening column density) are chosen. At the same time, the fits also suggest 
that their masses are not extremely large, generally < 150 - 300 M . Regrettably, however, 
our confidence that the formally best constrained of these (M81 XI) has a large mass must 
be tempered by the realization that the L/Le implied by even the highest acceptable mass 
for this object is ~ 0.85, a normalized luminosity so large as to call into question the physical 
standing of the model from which it was derived. 

It is also of interest that the three objects with the strongest lower bound on the inferred 
mass are exactly the three objects in our sample with the greatest luminosity, from ~ 3 - 
8 x 10 39 erg s _1 . Although M101 X2 and NGC 253 XI do not have inferred luminosities in 
Eddington units as high as that of M81 XI, they are not far below: ~ 0.3 - 1 for M101 X2 
and ~ 0.35 - 0.55 for NGC 253 XI. 



The remaining two (NGC 253 X3 and X4) have best-fit masses that are several tens 
of solar masses, but their data are also consistent with a mass both considerably lower 
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(< 10 M Q ) and considerably higher (> 100 M ). Circinus X2, as we have already discussed, 
is also in this category, but with an even wider range of uncertainty. All three of these have 
luminosities low enough, 3 - 10 x 10 38 erg s _1 , that they are at best borderline ULXs in any 
case. 

Along with the mass and inclination angle, our fits are also sensitive to the black hole 
rotation rate. Viewed in the mass-spin plane, \ 2 rises sharply in the case of the three best- 
constrained objects, M 81, M 101, and NGC 253 X-l when a/M < 0.8-0.9. There is a slight 
degeneracy between these two variables, in the sense that lower spin requires smaller mass, 
but it is quite weak. However, we do not regard these results as robust because of a potential 
systematic error: our assumption of a Novikov-Thorne surface brightness profile. 



If, as recent work on MHD stresses in accretion disks suggests ( IKrolik et al.l 120051 ) . 
there are significant stresses throughout the region of marginally stable (or even unstable) 
orbits, additional dissipation there is also likely. The result would be, for fixed luminosity, 
a shift to higher temp erature and a smaller effective radiating area in the emitted spectrum 



(lAgol fc Krolik 1 120001 ). Employing an estimate of the local dissipation rate derived from 



simulation data and a general relatiy i stic ra y-tracing code to find what portion of the radia- 



tion reaches infinity, iBeckwith et al.l (120081 ) have recently computed the effective "radiation 
edge" due to these effects. They found that when a/M < 0.9 and the inclination angle is 
such that cosz < 0.7, the characteristic radius for the radiation seen by distant observers can 
be displaced inward by factors of several relative to the prediction of the Novikov-Thorne 
model. 

The sense of the bias induced by fitting a spectrum on the basis of a Novikov-Thorne 
surface brightness profile if the radiation edge really is located closer to the black hole is 
to find a mass smaller and/or a spin larger than the actual one. A factor of 2 error in 
the characteristic radius in gravitational units (i.e., r/r g ) translates directly into a factor 
of 2 in mass because r g = GM/c 2 . In terms of spin, a factor of 2 error would mean that 
a nominal best-fit of a/M = 0.9 should be reinterpreted as actually indicating a/M = 0.4 
if the true radiation edge, like the radiation edge predicted by the Novikov-Thorne model, 
scales linearly with the radius of marginal stability. 

Given the several order of magnitude a priori uncertainty in our knowledge of the masses 
of ULXs, a factor of 2 systematic error in mass is, at this stage, not a major handicap. More- 
over, because we are especially concerned with determining whether ULX black hole masses 
can be as small as those found in Galactic black hole binaries, the fact that any correction 
would increase the inferred black hole mass means that the Novikov-Thorne estimate is a 
conservative one. 
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On the other hand, given the much more limited range of possible black hole rotation 
rates, a systematic error of order unity is a major concern in this context. In addition, the 
character of the spectral fitting acts to enhance the magnitude of the possible systematic 
error. Diminishing the spin in a model fit, by moving the radiation edge outward, also lowers 
the characteristic temperature for fixed luminosity and black hole mass. In order for the 
spectrum that results to fit the data, a larger orbital blue shift is required, fo rcing the model 



towar d greater inclination angle. Large inclination angles are exactly where iBeckwith et al. 



(120081 ) found the largest offsets between the radiation edge as predicted by the Novikov- 
Thorne model and the edge as inferred from simulation data. 



7. Summary 

After selecting the 23 known ULXs with the highest signal-to-noise and least absorbing 
column density, we chose a subsample of 6 objects in which the spectrum appeared to be 
purely that of a thermal disk, with no hint of any power- law component. Fitting their spectra 
with two different disk models (one based on the multi-color disk approximation, the other 
resting on detailed stellar atmosphere calculations), we found that in all but one case, the 
results were very similar: the masses yielding the best fits to the data lie between 20 and 
80 M Q . In the exception, the two models suggested very different masses, 340 M in the 
case of the atmosphere model, 13 M in the case of the multi-color disk model. 

More significantly, in 3 of the 6, the model fitting (that is, a search for those parameters 
for which a disk model can reproduce the observed spectrum to within the errors) formally 
excludes masses similar to those seen in Galactic black holes, ~ 10 M . The models clearly 
prefer rather larger masses, from several tens of solar masses to ~ 100 M Q , although in one 
case the value of L/Le preferred by the model-fitting is so large (~ 1) as to be physically 
inconsistent with a disk model. We also see formal upper bounds on the mass that are 
generally ~ 150 - 300 M , although systematic error having to do with the detailed surface 
brightness profile of relativistic disks may relax these bounds by a factor of 2 - 3. In the 
other three objects, the error bars are too broad to permit confident exclusion of either 
conventionally small masses or considerably larger ones. 

On this basis, we believe that these data provide significant new evidence that at least 
some ULXs have masses rather greater — ~ 30 - 100 M — than is found in ordinary 
Galactic binary black holes. 

We are grateful to Andy Ptak for extensive instruction in the proper treatment of X-ray 
spectral data. We also thank Shane Davis and Omer Blaes for many helpful discussions 
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about the calculation of disk atmospheres. We especially acknowledge Ivan Hubeny for 
construction and maintenance of the stellar atmosphere code TLUSTY and its disk version 
TLUSDISK. 

This work was partially supported by NASA ATP Grant NAG5-13228 and by NSF 
Grant AST-0507455. 
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Table 1. The Disk-Dominated ULX Sources 



Ob ID 


Galaxy 


D (Mpc) 


Object 


RA (h m s) 


Dec (o / //) 


Tot. Red. Cts 


Red. Ct. Rate (10~ 2 cts t" 1 ) 


0200980101 


M81 


3.6 1 


xl 


09 55 32.9 


+69 00 34.8 


6913, 6524, - 


9.12, 8.00, - 


0104260101 


M101 


7.4 2 


x2 


14 03 03.8 


+54 27 37 


874, 925, 1879 


3.22, 3.04, 9.55 


0152020101 


NGC 253 


3.73 3 


xl 


00 47 32.8 


-25 17 52.6 


5874, 6316, 16454 


8.78, 9.24, 28.46 








x3 


00 47 35.2 


-25 15 13.8 


2266, 2360, 5089 


3.33, 3.41, 8.72 








x4 


00 47 23.3 


-25 19 06.5 


801, 736, 1659 


1.16, 1.05, 2.82 


0111240101 


Circinus 


4.0 4 


x2 


14 12 54.2 


-65 22 55.3 


1391, -, 2703 


1.51, -, 4.07 



1 Frccdman ct al. i mil). 

2 iKelson et~aH [|l996h ; ljuvcevic fc Butcher] ||2006| ). 

3 Calculated from the distance modulus given in NED. Also refer to Karachcntscv ct al. | <2003ft ; lRekola et al.l 1120051) . 
4 lFreeman et alj fl977h ; llaria et alj feOOSh . 
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Table 2. Fitting Results 



Galaxy 


Object 




n H a 


M/IOOMq 


L/L e h 




Spin c 


Cos i d 


X 2 /dof 


Lobs/Lbol" 


M81 


xl 


1 


qn +0.l2 
yu -0.11 


u - '^-o.oe 


1.00-0.15 


0. 


997-o.oi 


n t;o+0.03 


386.3/397 


66.2/76.5 






1 


61 +o.n 

D1 -0.11 


4O+0.25 


o.82t° o: « 

0.75+H 


0. 


999_o.o5 


52+° 12 
u - oz -0.18 


392.6/397 




M101 


x2 


1 


18 -0.33 


n 57+ 1 - 21 

u -°'-0.27 


0. 


997-0.20 


n 5Q+° 15 

u - oa -0.38 


156.0/156 


41.5/46.1 









qfi +0.18 
yo -0.15 




0.38_ ; 26 





999-2.0 


n 4fi+°- 54 

u - 4D -0.37 


157.6/156 




NGC 253 


xl 


1 


fi o+0.12 
' bS -0.15 


72+ - 09 
u - ,z -0.22 


45+ - 11 

u '^°-0.10 





■997_„.i2 


34+ - 10 
u -° -0.04 


998.5/963 


29.2/31.5 






1 


cri+0.08 
ol -0.10 


0.73t° ;S 


23+ 08 
u - ZJ -0.05 


0. 


999-0.03 


n 90+0.09 

U - 2S -0.15 


999.3/963 






x3 


2 


qo+0.15 
• ac> -0.19 


0.50i°;^ 


u - zu -0.07 


0. 


997-0.20 


51 + ' 20 
u - 01 -0.17 


379.6/388 


10.1/11.0 






2 


7O+0.16 
la -0.13 


n 63 +113 


0Q+ - 09 
u - uy -0.06 


0. 


999_o.20 


n 42+0-28 
u '^ z -0.33 


380.5/388 






x4 





64 +0.21 


n 9 o+0.68 
u - zo -0.13 


14+ - 10 
u - i ^-0.08 





997_o.34 


46+ ' 27 
u -™-0.42 


143.9/140 


3.0/3.4 









45 +o.i8 

•^°-0.27 


n 35+ - 42 

u -°°-0.30 


n n5 +0 - 30 





999+»° 01 


n on+0.67 
U.JJ_ 25 


145.4/140 




Circinus 


x2 


5 


42 +0.27 
*-0.38 


3 40+ 12 - 70 


04 +0 - 60 
u - u ^-o.oi 


0. 


996+; : ° 01 


n n7+°- 68 

u - u '-0.07 


176.4/173 


7.3/8.6 






5 


45+0-46 


n 1 Q+3-28 


u - d4 -0.33 





501+?'' 498 


1 0O+ 000 


177.6/173 





a Column density in units of 10 21 cm -2 . 

b Luminosity normalized to the Eddington luminosity. 

c Dimcnsionlcss spin parameter of the black hole. 

d i - the inclination angle, i.e., the angle between the disk normal and the line of sight. 

e The ratio of the observed luminosity to the model integrated luminosity. Luminosity is in units of 10 38 erg s - 



Note. — 1) For each object, the first row shows the fitting results from BHSPEC and the second row with ". . ." 
shows the results from KERRBB; 2) All uncertainties were computed with A\ 2 = 2.706, equivalent to 90% confidence 
for a single parameter. 3) Where no upper uncertainty is shown, the data are consistent with all parameters up to the 
edge of the parameter space, i.e., for L/Le the upper bound is unity while for spin it is 0.997 for BHSPEC and 0.999 
for KERRBB. 
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channel energy (keV) 



channel energy (keV) 



Fig. 1. — Observed spectra and best-fit models (only BHSPEC shown). Panels from top to 
bottom: upper left — M81 XI, upper right — M101 XI, middle left — NGC 253 XI, middle 
right — NGC 253 X3, lower left — NGC 253 X4, lower right — Circinus X2. In each panel, 
the green, blue and red colors represent the data from MOS1, MOS2 and PN, respectively. 
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BHSPEC 



150 - 



KERRBB 



CO 

O 



100 - 



50 



M81 X1 M 1 01 X2 N253 X1 N253 X3 N253 X4 

OBJECTS 



Fig. 2. — The distribution of black hole masses in our sample (except Circinus X2). Masses 
derived using BHSPEC are shown by stars, those inferred on the basis of KERRBB are 
shown with open diamonds. All uncertainties were computed with A% 2 = 2.706, equivalent 
to 90% confidence for a single parameter. 
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Fig. 3. — Confidence level contours in the M-cosi plane. The location of the least x 2 is 
marked by a plus sign. A pair of dashed lines sh o ws th e extreme Kerr and Schwarzschild 
degeneracy relations suggested by ISun &: Malkanl (119891 ) . A dotted vertical line at M — 
20 Mm marks the greatest mass c onsistent with any of the known Galactic black hole binaries 
(jMcClintock &: Remillard II2004I ). The four contours show values of x 2 above the minimum 
by Ax 2 = 2.30 (red), 6.14 (green), 9.21 (blue), 13.82 (magenta). These correspond to the 
68.3% (1-0-), 95.4% (2-a), 99.0% and 99.9% confidence levels, respectively. Short ticks along 
each contour show the "downhill" direction in x 2 ■ 



